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The cosmological evolution of neutral hydrogen is an efficient way of tracing struc- 
ture formation with redshift. It indicates the rate of evolution of gas into stars 
and hence the gas consumption and rate star formation history of the Universe. In 
measuring HI, quasar absorbers have proven to be an ideal tool and we use obser- 
vations from a recent survey for high-redshift quasar absorption systems together 
with data gathered from the literature to measure the cosmological comoving mass 
density of neutral gas. This paper assumes f2jy/=0.3, £1/^=0.7 and h=0.65. 

1 Quasar Absorption Systems 

Quasar absorbers are systems along the line-of-sight between an observer and a 
quasar, and thus are observed in absorption in the spectra of the background 
source. They are classified according to their HI column density: Damped 
Lyman-a systems (hereafter DLA) have N(HI) > 10 20 ' 3 atom cm~ 2 , Ly- 
man Limit Systems (hereafter LLS) have N(HI) > 10 172 atom cm -2 and 
the lyman-a forest is composed of systems with smaller column densities. 
The work presented here is based on observations from a recent high-z sur- 
vey (Peroux et al. 2001) together with a compilation of absorption systems 
from the literature (Storrie-Lombardi and Wolfe 2000). This included 82 
quasars with z > 4 and resulted in a sufficiently large sample of quasar ab- 
sorption systems to enable the study of their properties with redshift (z). Fig 
1 shows the cumulative number of quasar absorbers per unit dX(z), for two z 
ranges, absorbers with z < 3.5 and z > 3.5, where dX is the distance interval: 
X(z) = + zf [(l + z) 2 (l + zn M )-z(2 + z)n A ]~ 1/2 dz. The column 

density of LLS cannot be derived directly from their equivalent widths as 
they correspond to the non-linear part of the curve of growth. Nevertheless, 
the expected number of LLS can be calculated by fitting N a (l + z) 1 to the 
observed number of systems per unit z and integrating over the z-path sur- 
veyed: LLS eX p ec t e d = J27=i I z ""^ + z )~' dz . The observations presented 
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Figure 1. Cumulative number of absorbers per unit distance interval, dX, for two different 
2 ranges. The "steps" correspond to the observed number of DLAs and the data points at 
log N(H I) = 17.2 are the expected number of LLS derived from the observed number of 
LLS per unit z. T-functions are used for the fits to the data. 



in Fig 1 show that there are fewer systems at the highest column density for 
z > 3.5 compared to the lower redshift range, which indicates that we are 
probing the epoch of assembly of high column density systems from smaller 
column density units. 

2 Cosmological Evolution of Neutral Hydrogen 

The neutral gas mass is calculated by integrating the quasar absorbers column 
density distribution, f(N, z): Q. HI {z) = Nf(N, z)dN . We used 

a T- function to fit the data (Fig 1) in order to take into account the HI 
included below the formal DLA definition. Indeed, DLAs were previously 
thought to contain the majority of the total HI but our new study shows 
that 14% at z < 3.5 and up to 40% at z > 3.5 of the neutral gas is actually 
contained in smaller column density systems (10 17 2 < N(HI) < 10 20 3 ). The 
total HI decreases with decreasing z thus matching the pattern expected from 
gas depletion due to star formation. Details on the method used and more 
discussion on the results presented here will be available in a subsequent paper. 
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Figure 2. Cosmological evolution of HI. The circles show the neutral gas contained in DLAs 
(filled circles are this work, open circles are from Rao and Turnshek 2000 who used HST 
spectra of quasars with know Mgll systems to look for low-z DLAs. Their results are 
surprisingly high but might be biased by small number statistics or the effect of gravitational 
lensing.) Vertical error bars correspond to lcr uncertainties and the horizontal error bars 
indicate bin sizes. The stars are the total HI including a correction for the neutral gas 
contained in systems with column densities below the DLA formal definition. The triangle 
at z = is the local HI mass measured by Rao & Briggs 1993. The squares, Qpnp 
and f2<3o (Fukugita, Hogan & Peebles 1998 and Gnedin & Ostriker 1992 respectively) are 
tlbaryons m local galaxies. Semi-analytical models from Somerville, Primack & Faber 2000 
are overplotted. 
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